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Q . Abstract. The low Mach number limit for the full compressible magnetohy- 

drodynamic equations with general initial data is rigorously justified in the 
whole space . The uniform estimates of the solutions in Sobolev space are 
^ obtained on a time interval independent of the Mach number. The limits are 

proved by using a theorem of G. Metiver & S. Schochet that established the 
decay of energy of the acoustic equations. 
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1. Introduction 

In this paper we study the low Mach number limit of local smooth solutions 

to the following full compressible magnetohydrodynamic (MHD) equations with 
general initial data in the whole space M.^ (sec [20,27,36,42]): 

dtp + divipu) ^ 0, (1.1) 

9t(/9u)+div(pu®u) + VP= -^-(curlH) x H + div4'(u), (1.2) 

An 

StH - curl (u X H) = -curl (i/ curl H), divH = 0, (1.3) 

dt£ + div (u(f ' + P)) = — div((u X H) X H) 
47r 
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+ div(^ — Hx (curlH) + u«'(u) + KV6lj. (1.4) 

Here the unknowns p, u = (ui, M2, ua) G M'^ , H = i?2, ^^3) G K'^, and 9 denote 
the density, velocity, magnetic field, and temperature, respectively; ^'(u) is the 
viscous stress tensor given by 

*(u) = 2/xro(u) + Adivu I3 

with ©(u) = (Vu + Vu~'')/2, I3 the 3x3 identity matrix, and Vu""" the transpose 
of the matrix Vu; S is the total energy given by £" = f + |Hp/(87r) and S' = 
p(e+|up/2) with e being the internal energy, p|up/2 the kinetic energy, and 
|Hp/(87r) the magnetic energy. The viscosity coefficients A and ^ of the flow 
satisfy /i > and 2/j, + 3A > 0. The parameter v > \s the magnetic diffusion 
coefficient of the magnetic field and k > the heat conductivity. For simplicity, we 
assume that /z. A, v and n are constants. The equations of state P = P{p, 6) and 
e = e(/5, 0) relate the pressure P and the internal energy e to the density p and the 
temperature of the flow. 

Multiplying (1.2) by u and (1.3) by H/(47r) and summing over, one finds that 

= div^t • u + — (curl H) x H • u + —curl (u x H) • H 

Att Att 

- —curl (curl H) -H. (1.5) 

Due to the identities 

div(H X (curlH)) = IcurlHp - curl(curlH) • H, 

div((u X H) X H) = (curlH) x H • u + curl (u x H) • H, (1.6) 

one can subtract (1.5) from (1.4) to rewrite the energy equation (1.4) in terms of 
the internal energy as 

dtipe) + dhilpvLe) + (divu)P = — IcurlHP + *(u) : Vu + kA6', (1.7) 

Att 

where ^'(u) : Vu denotes the scalar product of two matrices: 

*(u) : Vu = ^ ^ ("1^ + + A|divu|2 = 2/i|ID)(u)|2 + A|trD(u)r 

To establish the low Mach number limit for the system (1.1)-(1.3) and (1.7), in 
this paper we shall focus on the ionized fluids obeying the perfect gas relations 

P = d\pd, e = eve, (1.8) 

where the parameters 91 > and cy > are the gas constant and the heat capacity 
at constant volume, respectively, which will be assumed to be one for simplicity of 
the presentation. We also ignore the coefficient l/(47r) in the magnetic field. 

Let e be the Mach number, which is a dimensionless number. Consider the 
system (1.1)-(1.3), (1.7) in the physical regime: 

P-Po + O(e), u-O(e), H - 0(e), V0 - 0(1), 

where Pg > is a certain given constant which is normalized to be Pq = 1- Thus we 
consider the case when the pressure P is a small perturbation of the given state 1 
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while the temperature 9 has a finite variation. As in [2] , we introduce the fohowing 
transformation to ensure the positivity of P and 9 

P{x,t) = e^P'^'^'^'l 9{x,t) = e^'^^'''\ (1.9) 

where a longer time scale t = r/e (still denote r by t later for simplicity) is intro- 
duced in order to seize the evolution of the fluctuations. Note that (1.8) and (1.9) 
imply that p{x,t) = e^p'(=^.<^*)-s'(^^^*) since = = 1. Set 

H(x,i) = eH"(a;,ei), u{x,t) = evL\x,et), (1.10) 

and 

^ = e//*^, A = eA^, v = ei/*^, k = en^. (l-H) 

Under these changes of variables and coefficients, the system, (1.1)~(1.3), (1.7) 
with (1.8), takes the following equivalent form: 

dtp' + (u^ • V)/ + idiv(2u^ - K^e-^P'+'^'ve^) 

= ee-'P'[j/'|curlH^|2 + ^'(u^) : Vu^] + n'e-'P'+^'Wp' ■ W9\ (1.12) 

e-^'[atU^ + (u^-V)u^] + — = e-^P'[(curlH^) x H"^ + div*^(u^)], (1.13) 

dtiV - curl (u' X H'=) - i^'^AH' = 0, divH' = 0, (1.14) 
819" + (u' • \/)9' + divu' 

= e2g-£p'[j,e|curlH'^|2 + ^'<=(u<=) : Vu<=] + K<=e-'^P'div(e^' V6i<^), (1.15) 

where ^'^(u^) 2^^D(u') + AMivu"^ I3, and the identity curl (curlH') = VdivH' - 
AH*^ and the constraint divH' = are used. 

We shall study the limit as e — > of solutions to (1.12)-(1.15). Formally, as 
e goes to zero, if the sequence (u'^, H', 6''^) converges strongly to a limit (w, B,i?) 
in some sense, and (/i*^, A"^, i^*^, k*^) converges to a constant vector (/l, A,P, k), then 
taking the limit to (1.12)-(1.15), we have 

div(2w- Ke^'Vz?) = 0, (1.16) 

e^'^[dtyv + {-w ■ + Vtt = (curlB) x B + div$(w), (1.17) 

dtB - curl (vi^ X B) - t?AB = 0, divB = 0, (1.18) 

dtd + (w • V)d + divw = Kdiv(e''V?9), (1.19) 
with some function tt, where $(w) is defined by 

$(w) = 2^D(w) + Adivwig. (1.20) 

The purpose of this paper is to establish the above limit process rigorously. For 
this purpose, we supplement the system (1.12)-(1.15) with the following initial 
conditions 

(p^u^H^r)|,=o = (pf„(a;),uf„(x),Hf„(x),0f„(a;)), x&R\ (1.21) 

For simplicity of presentation, assume that /i'^ = /I > 0, z^'^ = > 0, k*^ = k > 0, 
and A*^ = A. The general case /i*^ — >■ /2 > 0, z^*^ -> i> > 0, — )• k > and A' — ?► A as 
e can be treated by modifying slightly the arguments presented here. 
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As in [2], we will use the notation ||u||/f^<^ := + ^yllwH/f" for any cr e M 

and > 0. For each e > 0, t > and s > 0, we will also use the following norm: 

\\ip^,u\ll\0^-e){t)\Ue 

:= sup {||(p^u^H^)(r)|U.-+|l(ep^eu^eH^0^-0)(r)|!^.+.} 

re[0,t] 

+ {^*[|lV(p^u^H^)||?,.+||V(eu^^i^^n|i^.+.](r)dr}'^^ (1.22) 

Then, the main result of this paper reads as follows. 

Theorem 1.1. Let s > 4. Assume that the initial data (pf^^^uf^jUf^jOf^) satisfy 

\\{pl,ul,Ul,ei^-m)\\s.c<Lo (1.23) 

for all e € (0,1] and two given positive constants 9 and Lq. Then there exist positive 
constants Tq and eo < 1, depending only on Lq and 9 , such that the Cauchy problem 
(1.12)-(1.15), (1.21) has a unique solution (p*^, u*^, H*^, 0*^) satisfying 

\\{p\n\YL\e' -e){t)\\,,,<L, Vte [0,ro], Vee (0,eo], (1.24) 

where L depends only on Lq, 8 and Tq. Moreover, assume further that the initial 
data satisfy the following conditions 

mx)~9\<No\x\-^-^, \S/e'o{x)\<No\x\~'''i, VeG(0,l], (1.25) 

(K„, curl (e-^'-ufJ, Hf„, 9^) ^ (0, wo, Bq, ^o) m H-\M.^) as e ^ 0, (1.26) 

where Nq and C are fixed positive constants. Then the solution sequence ij)^ , u*^, H*^, 
9^) converges weakly in L°°{0,To; H^R^)) and strongly in 1/^(0, Tq; H^^^{R^)) for 
all < S2 < s to the limit (0,w, B,i?), where (w, B,i9) satisfies the system (1.16)- 
(1.19) with initial data (w, B,i?)|(=o = ("Wq, Bq, i^o)- 

We now give some comments on the proof of Theorem 1.1. The key point in 
the proof is to establish the uniform estimates in Sobolev norms for the acoustic 
components of solutions, which are propagated by wave equations whose coefficients 
are functions of the temperature. Our main strategy is to bound the norm of 
(Vp*^, divu') in terms of the norm of {edt){p'^, u*^, H*^) and that of {ep'^, eu'^, eH'^, 9'^) 
through the density and the momentum equations. This approach is motivated by 
the previous works due to Alazard in [1,2], and Levermore, Sun and Trivisa [37]. 
It should be pointed out that the analysis for (1.12)-(1.15) is complicated and 
difficult due to the strong coupling of the hydrodynamic motion and the magnetic 
fields. Moreover, it is observed that the terms (curlH^) x in the momentum 
equations, curl (u*^ x H^) in the magnetic field equation, and |V x H'^p in the 
temperature equation change basically the structure of the system. More efforts 
should be payed on the estimates involving these terms, in particular on the higher 
order spatial derivatives. We shall exploit the special structure of the system to 
obtained the tamed estimate on higher order derivatives, so that we can enclose 
our arguments on the uniform boundedness of the solutions. Once the uniform 
bounds of the solutions are obtained, the convergence result in Theorem 1.1 can 
be proved by applying the compactness arguments and the dispersive estimates on 
the acoustic wave equations in the whole space developed in [39] . 

Remark 1.1. The positivity of the coefficients /x, v and k plays an fundamental role 
in the proof of Theorem 1.1. The arguments given in this paper can not be applied 
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to the case when one of them disappears. We shall discuss this situation in another 
forthcoming paper. 

We point out that the low Mach number limit is an interesting topic in fluid 
dynamics and applied mathematics. Now we briefly review some related results 
on the Euler, Navier-Stokes and MHD equations. In [47], Schochet obtained the 
convergence of the non-isentropic compressible Euler equations to the incompress- 
ible non-isentropic Euler equations in a bounded domain for local smooth solutions 
and well-prepared initial data. As mentioned above, in [39] Metivier and Schochet 
proved rigorously the incompressible limit of the compressible non-isentropic Euler 
equations in the whole space with general initial data, see also [1,2,37] for further 
extensions. In [40] Metivier and Schochet showed the incompressible limit of the 
one-dimcnsional non-isentropic Euler equations in a periodic domain with general 
data. For compressible heat-conducting flows, Hagstrom and Lorenz established 
in [18] the low Mach number limit under the assumption that the variations den- 
sity and temperature are small. In the case of without heat conductivity, Kim and 
Lee [32] investigated the incompressible limit to the non-isentropic Navier-Stokes 
equations in a periodic domain with well-prepared data, while Jiang and Ou [31] 
investigated the incompressible limit in three-dimensional bounded domains, also 
for well-prepared data. The justification of the low Mach number limit to non- 
isentropic Euler or Navier-Stokes equations with general initial data in bounded 
domains or multi-dimensional periodic domains is still open. We refer the inter- 
ested reader to [6] for formal computations for viscous polytropic gases, and to [5,40] 
for the study on the acoustic waves of the non-isentropic Euler equations in periodic 
domains. Compared with the non-isentropic case, the description of the propaga- 
tion of oscillations in the isentropic case is simpler and there are many articles on 
this topic in the literature, see, for example, Ukai [49], Asano [3], Desjardins and 
Grenier [11] in the whole space case; Isozaki [25,26] in the case of exterior domains; 
Iguchi [24] in the half space case; Schochet [46] and Gallagher [16] in the case of 
periodic domains; and Lions and Masmoudi [43], and Desjardins, et al. [12] in the 
case of bounded domains. 

For the compressible isentropic MHD equations, the justification of the low Mach 
number limit has been established in several aspects. In [33] Klainerman and Majda 
studied the low Mach number limit to the compressible isentropic MHD equations in 
the spatially periodic case with well-prepared initial data. Recently, the low Mach 
number limit to the compressible isentropic viscous (including both viscosity and 
magnetic diffusivity) MHD equations with general data was studied in [23,28,29]. 
In [23] Hu and Wang obtained the convergence of weak solutions to the compressible 
viscous MHD equations in bounded domains, periodic domains and the whole space. 
In [28] Jiang, Ju and Li employed the modulated energy method to verify the limit 
of weak solutions of the compressible MHD equations in the torus to the strong 
solution of the incompressible viscous or partially viscous MHD equations (zero 
shear viscosity but with magnetic diffusion), while in [29] the convergence of weak 
solutions of the viscous compressible MHD equations to the strong solution of the 
ideal incompressible MHD equations in the whole space was established by using 
the dispersion property of the wave equation, as both shear viscosity and magnetic 
diffusion coefficients go to zero. For the full compressible MHD equations, the 
incompressible limit in the framework of the so-called variational solutions was 
established in [34,35,41]. Recently, the low Mach number limit for the ideal and full 
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non-iseiitropic MHD equations with small entropy or temperature variations was 
justified rigourously in [30]. Besides the references mentioned above, the interested 
reader can refer to the monograph [14] and the survey papers [9, 44, 48] for more 
related results on the low Mach number limit to fluid models. 

We also mention that there arc a lot of articles in the literatures on the other 
topics related to the compressible MHD equations due to theirs physical importance, 
complexity, rich phenomena, and mathematical challenges, sec, for example, [4, 7, 
8, 10, 13, 15, 19-22, 38, 42, 51] and the references cited therein. 

This paper is arranged as follows. In Section 2, we describe some notations, recall 
basic facts and present commutators estimates. In Section 3 we first establish a 
priori estimates on (H'^, 9'^), {ep'^, eu'^, eH^, 9'^) and on {p'^, u'^). Then, with the help 
of these estimates we establish the uniform boundeness of the solutions and prove 
the existence part of Theorem 1.1. Finally, in Section 4 we study the decay of 
the local energy to the acoustic wave equations and prove the convergence part of 
Theorem 1.1. 

2. Preliminary 

In this section, we give some notations and recall basic facts which will be used 
frequently throughout the paper. We also present some commutators estimates 
introduced in [37] and state the results on local solutions to the Cauchy problem 
(1.12H1.15), (1.21). 

Wc denote (•, •) the standard inner product in L^(]R'^) with norm (/, /) — 
11/11^2 and H'' the standard Sobolev space W'^''^ with norm || • The notation 
II {Ai, . . . , y4.fc)||^2 means the summation of || Ai||^2, i = 1, . . . , fc, and it also applies 
to other norms. For the multi-index a = (ai,a2,a3), we denote 9" = d^^d^^d^^ 
and |a| — \ai\ + \a2\ + [asl. Wc will omit the spatial domain M."^ in integrals for 
convenience. We use k > {i £ N) to denote given constants. We also use the 
symbol K or Co to denote generic positive constants, and C(-) to denote a smooth 
function which may vary from line to line. 

For a scalar function /, vector functions a and b, we have the following basic 
vector identities: 



div(/a) 


^ /diva + V/ • a. 


(2.1) 


curl (/a) 


= / • curia — V/ X a. 


(2.2) 


div(a X b) 


= b • curl a — a • curl b. 


(2.3) 


curl (a X b) 


= (b • V)a - (a • V)b + a(divb) - b(diva). 


(2.4) 


V(a-b) 


= (a • V)b + (b • V)a + a x (curlb) + b x (curia). 


(2.5) 



For simplicity of presentation, we will use 

Definition 2.1 ( [2]). For s e and g e H-''+'^{R^), wc define the weighted 
Sobolev norm || • as 

llffllffl+i = llffllff-' +e||V.9||H^-. 
Below wc recall some results on commutators estimates. 

Lemma 2.1. Let a = (ai,a2,a3) be a multi-index such that \a\ ~ k. Then, for 
any cr > 0, there exists a positive con stant Cq, such that for all f,g £ H''+'^{R^), 

||[/,a"].9||H^ <Co(!|/|ki,.cH,9||^„+.-i + ||/||ff.+.||.9||L^0. (2.6) 
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Lemma 2.2. Let s > 5/2. Then there exists a positive constant Co, such that for 
all e e (0, 1], T > and multi-index (3 = ^2, /?3) satisfying < |/?| < s — 1, and 
any f,g e C°°([0, T], iJ^(K3)), it holds that 

\\[f,d^{edt)]gh2 <eCoi\\f\\H^^-4dtg\\H^-^ + \\dtf\\H^-4g\\H^,-i). (2.7) 

Since the system (1.1)-(1.3), (1.7), (1.8) is hyperbolic-parabolic, so the classical 
result of Vol'pert and Khudiaev [50] implies that 

Proposition 2.3. Let s > 4. Assume that the initial data (po, Uq, Hp, ^o) satisfy 

\\{po - p, uo,Ho,6'o -d)\\H= < Co 

for some positive constants p, 9 and Co- Then there exists a T > 0, such that the 
system (1.1)-(1.3), (1-7), and (1.8) with these initial data has a unique classical 
solution {p,u,H,e) enjoying pe C{[0,f],H''{R^)), (u, H, 6*) e C{[0,f], H'{R^)) n 
L^{0,f;H''+^R^)), and 

sup II (p,u,H, 0)111,.+ r{Ai||D(u)||l,. +A||divu||l,. 

+ j/||VH||2,. +k||V0||^3} (r)dr < iCl 

It follows from Proposition 2.3 and the transforms (1.9), (1-10) that there ex- 
ists a > 0, depending on e and Lo, such that for each fixed e and any ini- 
tial data (1.21) satisfying (1.23), the Cauchy problem (1.12)-(1.15), (1.21) has a 
unique solution (p^u^H^6l^) satisfying {p^ ,u' ,W ,0") e C{[0,T,), H'iR^)) and 
(u%H%6'^) n i2(o,r,;iJ'*+i(R3)). Moreover, let T* be the maximal time of exis- 
tence of such smooth solution, then if T* is finite, one has 

hmsup {||(p^u^H^0^)(^)||H.l.=o + ||(u^H^0^)(^)||H'-,=o} = oo. 

Therefore, the existence part of Theorem 1.1 is a consequence of the above 
assertion and the following a priori estimates which can be shown in the same 
manner as in [39]. 

Proposition 2.4. For any given s > 4 and fixed e > 0, let (p*^, u"^, H*^, 6*^) be the 
classical solution to the Cauchy problem (1.12)-(1.15) and (1.21). Denote 

O(T):=||(g^u^H^0^-^)(T)||,,„ (2.8) 

Oo := ||(Kn,uL,Hf„,0fJ||«. + ||e(Kn,uf„,Hf„,0f„ - 0)\\h;+-- (2-9) 

Then there exist positive constants To and eo < 1, and an increasing positive func- 
tion C(-), such that for all T S [0,To] and e G (0,60]; 

0{T) < C{Oo)exp{{VT + e)C{0{T))}. 

3. Uniform estimates 

In this section we shall establish the uniform bounds of the solutions to the 
Cauchy problem (1.12)-(1.15) and (1-21) by modifying the approaches developed 
in [2,37,39] and making careful use of the special structure of the system (1.12)- 
(1.15). In the rest of this section, we will drop the superscripts e of the variables 
in the Cauchy problem and denote 

4'(u) = 2/2D(u) + Adivu I3. 
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Recall that it has been assumed that = p, > 0, v'^ = D > 0, k'^ = R > 0, and 
A*^ = A independent of e. 

3.1. i/* estimates on (H,6') and {ep,eu). To prove the Proposition 2.4, we first 
give some estimates derived directly from the system (1.12)-(1.15). Denote 

Q : = II (p, u, H, - 0)\\h' + II {ep, eu, eH, 9 - 0)||^.+2, (3.1) 

S : = ||(Vu,Vp,VH)||ff= + ||V(eu,eH,0)||^.+2. (3.2) 

One has 

Lemma 3.1. There exists an increasing function C(-) such that 

sup {||(H,0)||h. + ||eH||^.+.} + ||V(H,0)|U.(o,T;H») + || V(6H)||^.(o 
te[o,T] 

< C'(C'o)cxp{TC(0)}. (3.3) 

Proof. For any multi-index a = (ai, 02,0:3) satisfying |q;| < s, let = c)"H. 
Then 

atH„ + (u • V)H„ - i>AH„ -[9", u] • VH - ^"(Hdivu) + 9"((H • V)u). 

Taking inner product of the above equations with and integrating by parts, we 
have 

i^||H„||i. + i>|| VH^IIi. = - ((u • V)H„, H„) - ([a", u] • VH, H„) 

- (9"(Hdivu), H„) + (a"((H • V)u), H„). 

An integration by parts gives 

-((u- V)H„,H„) = i J divu|H„|2dx<C(0)||H„||i.. 

It follows from the commutator inequality (2.6) that 

||[9",u] • VHIU2 < Co(||u||h.i,oo||VH||h-^-i + ||u||h. || VH|Uc.) < C{0). 
By Sobolev's inequality, one gets 

~(a"(Hdivu),H„) + (a"(H.Vu),H„) <Co||H||2,.||u||h=+i <C(0)5. 
Thus, we conclude that 

||H||ff. + i?||VH|U2(o,T;ff=) < C(Oo) + C{0)T + C{0) [ S{T)dT 

Jo 

< C(Oo) + C{0)T + C{0)Vt 
<C(Oo) + C(0)Vt 
<C(Oo)e^^('^). 

Now denote H = eH and Hq = 9"(eH) for jaj = s + 1. Then, Hq satisfies 

dtiia + (u • V)H„ - j>H„ = -e[d",u] ■ VH - ea"(Hdivu) + e9"((H • V)u). 

(3.4) 

The commutator inequality (2.6) implies that 

II -e[9",u] • VHIU2 < Co(||u||„,i,oo||VH|U. + ||eu||^.+i||VH||i^) < C{0), 
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while an integration by parts and Sobolev's inequality lead to 
-(ea"(Hdivu),H„) + (ea"((H. V)u),H„) < ^||VH„||2, + Co||H|||,.||eu||2,.+i 

<^||VH„|ii.+C(0). 

Hence, we obtain 

||eH||^=+i + j?|lV(eH)|U.(o,T;/f»+i) < C(Oo) + C{0)T. 

Similarly, the terms e^||Hj|/fs+2 and e^|| VH||i2(Q j-.^s+i) can be estimated. 
Next, we estimate 9. Using Sobolev's inequality, one finds that 

\\d''{e^e-'P[D\cuTlU\'^ + l'(u) : Vu])||i2 

< Co||ep|U.(||eH||H.+i + \\eu\\H^) < C{0). 
Employing arguments similar to those used for H, we can obtain 
II^IIh^ + ||V0|U2(o.T;ff») <C{Oo)exp{TC{0)}. 
Thus, the lemma is proved. □ 

Lemma 3.2. There exists an increasing function C(-) such that 

sup \\{ep,eu)\\Hs + ||V(eu)|U2(o,T;H=) < C(Oo) cxp{VTC(0)}. (3.5) 

Proof. Let p = ep, and pa = 9"(ep) for any multi-index a = (ai,a2,a3) satisfying 
|q;| < s. Then 

dtpa + (u • W)pa = - [9", u] • (Vp) - 9"[div(2u - Ka{ep)b{e)We)] 
+ a"{a(ep)[i/|curl (eH)p + *(eu) : V(eu)]} 
+ Kd"{a{ep)b{e)V{ep) ■ V9} 

+ h2 + h3 + hi, (3.6) 
where, for simplicity of presentation, we have set 

a{ep) e'^P, b{0) e^. 
It is easy to see that the energy estimate for (3.6) gives 

= -((u ■ V)25„,p„) + {hi + h2 + h3 + hi,Pa), (3.7) 

where we have to estimate each term on the right-hand side of (3.7). First, an 
integration by parts yields 

-((U • V)Pa,Pa) divulp^pdx < C{0)\\Pa\\l., 

while the commutator inequality leads to 

\\hi\\ < Co(||u||wi-||Vj5||h-i + \\u\\h4^P\\l^) < C{0). 
Consequently, 

{hl,Pa) < IIPalUHI/iilU^ < CiO). 
From Sobolev's inequality one gets 

||/^2|| < Co||u||h.+i + \\e\\Hs\\ep\\H^\\0\\H^^2 < C{S){1 + C{0)), 
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whence, 

{h2,Pc.) <C(0)C{S). 
Similarly, one can prove that 

((/i3 + M,Pa> <cio). 

Hence, we conclude that 

\\ep\\H^ <C{0„)cxp{VTC{0)}. 

In a similar way, we can obtain estimates on u. Thus the proof of the lemma is 
completed. □ 

Next, we control the term j|(u,p)||j:/s. The idea is to bound the norm of (divu, 
Vp) in terms of the suitable norm of (eu, ep, eH, 9) and e(9tu, dtp) by making use 
of the structure of the system. To this end, we first estimate |j (eu, ep, 6')||^.,+i . 

3.2. estimates on (eu, ep, eH, 0). Following [2], we set 

{p, u, H, 9) := {ep - 9, eu, eH, 9 - 9). 

A straightforward calculation results in that (p, u, H, 9) solves the following system: 

dtp + (u • V)p + idivu = 0, (3.8) 

h(-9)\dtu + (u • V)u] + i(V]5 + V^) 

= a(ep)[(curlH) x H + div*(u)], (3.9) 
atH + u-VH + Hdivu-H-Vu-PAH = 0, divH = 0, (3.10) 

dt9 + (u • V)^ + idivu = ea(ep)[z?curlH : curlH + ea(ep)^'(u) : Vu] 
e 

+ Ka(ep)div(&(6')V^). (3.11) 

One has 

Lemma 3.3. Let s > 4 and (p, u, H, 6') he a solution to the problem (1.12)-(1.15), 
(1.21) on [0,Ti]. Then there exists a constant li > 0, such that for any e € (0, 1] 
and t Cz [0, T], r = min{ri, 1}, it holds that 

sup \\ieq,eu,9-9)iT)\\Hs+i+h{ [ || V(eu, 0)||^,+i (T)dT| 

relO.T] '-Jo ^ 

<C{Oo)cxp{VtC{0{T))}, (3.12) 
where 0{T) and Oq are defined by (2.8) and (2.9), respectively. 
Proof. Let a = (ai, a2, 0^3) be a multi-index such that |a| = s + 1. Set 
(p„, u„. He., 9c.) ■■= (a"(ep - 9), a"(eu), a"(eH), 9"(0 - 9)) . 
Then, H^ satisfies (3.4) and (Pq,Uq,6'q) solves 

dtPa + {u-'V)pa + ^divua = gi, (3.13) 

b{-9)[dtUa + (u • V)u„] + i(Vp„ + V^„) 

= a(ep)(curlH) x H^ + a(ep)div*(u„) + 52, (3.14) 
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dtOa + (u • V)^Q, H — divUa ~ ea(ep)[S'CurlH : curlHa + ^'(u) : VUq] 

e 

+ Sa(ep)div(6(6l)V0„)+g3 (3.15) 

with initial data 

(pa,Ua,Hc<,0„)|t=o := (9"(epin(a;) - 6'in(a;)),9"(euin(a:;)), 

9"Hi„(x)),a"(0i„(x)-0")), (3.16) 

where 

51 :=-[a",u].V(ep-0), (3.17) 

52 [a",&(-0)]9t(eu) - [9", &(-0)u] • V(eu) 

+ [9", a(ep)curl (eH)] x H + [9", a(ep)]div«'(eu), (3.18) 

53 := - [9", u] • + [9", a(ep)curl (eH)] : curl (eH) 

+ e[a",a(ep)^'(u)] : V(eu) + (a(ep)div(6(6')V6')) 

- Ka(ep)div(6(6')V^a). (3.19) 

It follows from Proposition 2.3 and the positivity of a(-) and h{-) that a(-) and 
fe(-) are bounded away from uniformly with respect to e, i.e. 

a{ep) > a > 0, b{-e) > 6 > 0. (3.20) 

The standard L^-energy estimates for (3.13), (3.14) and (3.15) yield that 

1 d 



:(l|Pa|li2 + (6(-6')Ua,U„) + \\Ua\\L2) 



I I 7 I „ II T n -1— llll F/lll„ ll„ J -I- I I F / „ J ^ 

2dt 

< i(6t(0)u„,Ua) - ((u- V)j3„,p„) - (6(-(?)(u- V)u„,u„) - ((u • V)^„, 



2 

+ (a(ep)(curlH) x Hq,Uq) + (a(ep)divvI'(uQ), u„) 

(ea(ep)[z?curlH : curlH^ + «'(u) : VuQ],6'a) + (Ka(ep)div(6(6')V6'a), ^a) 
+ {gi,Pa) + {92,Ua) + {93,L)- (3.21) 
It follows from equation (1.15) and the definition of Q and S that 
WhrniL^ < \\b{9)\\H40t\\H^ < C{Q){1+S). 

Therefore, 

^{btie)u^,Ua) <C{Q){1+S). 

On the other hand, it is easy to see that 

-((u • V)p„,pa) - {b{-e){u ■ V)u„, Ua) - ((u • V)^„, < CiQ) 

and 

(a(ep)(curlH)xH,,u„) <C(Q). 
A partial integration shows that 

~(a(epo)div^'(u),u) J /^a(epo)[|Vu„p + (^ + A)|divu„|2]da; 

+ ^{{\/a{ep) ■ V)Uq,Uq) 

+ + A)((Va(ep)divuQ, UQ,)dx 
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= rfl+fi2 + rf3. (3.22) 

Thanks to the assumption that p, > and 2/2 + 3A > 0, there exists a positive 
constant , such that 

di>aC J iVuafdx, (3.23) 

while Cauchy-Schwarz's inequahty imphes 

\d2\ + \d3\<C{Q)S. (3.24) 

Similarly, one can obtain 

-{Ka{ep)div{b{0)V§a),9a) > Safe||V0„|||2 - C{Q)S. (3.25) 

Easily, one has 

|(ea(ep)[PcurlH : curlH„ + 1'(u) : Vuc.],^a)| < C(Q)(l+5). 

It remains to estimate {gi,Pa), (32, Mq) and (.93,^0) in (3.21). First, an applica- 
tion of Holder's inequality gives 

\{pa,gi)\ < CoWPah^hlh'^, 
where can be bounded, by using (2.6), as follows 

llffilU^ =||[a",u].V(ep-0)|U2 

<Co(i|u||wi,oo||V(ep-0)|U. + j|u||a.+i|lV(ep-0)|U.^.). 
It follows from the definition of Q and Sobolev's inequalities that 

\\y{ep,e)\\H^ < Q, ||V(ep,0)||L=o < Q. 
Therefore, we obtain HpiHl^ < C(Q)(1 +5), and 

\{Pa,gi)\<C{Q){l+S). (3.26) 
Next, we turn to the term |(uQ,, g2)|- Due to the equation (1.13), one has 
~[d",b{~0)]dt{eu) =[9^6(-^^)]((u• V)(eu)) + [9", 6(-0)] (&(0)Vp) 

- [d°',b{~d)]{b{e)a{ep){cm\U x (eH))) 

- [5", b{~e)] (6(e')a(ep)div*(eu)) . (3.27) 
The inequality (2.6) implies that 

|(u„,[a",&(-0)]((u.V)(eu)))| 
<Co\\ujL2\\[d'',b{-e)]{u-V{eu))\\L2 

< CiQ){\\bi-9)\\w^.^\\u ■ V(eu)||H., + ||6(-(?)||h.+i ||u • V(eu)||ioo 

< C{Q), 

and 

\{u^,[d^,b{-e)]{b{e)\/p))\ 
<Co\\ujL4[d",b{~e)]{b{e)Vp)\\L2 

< c{Q){\\b{-e)\\nn^o.\\b{9)Wp\\H-, + \\b{-e)\\H^^4b{e)Wp\\L^ 

<C(Q)(1 + 5). 
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The third term on the right-hand side of (3.27) can be treated in a similar manner, 
and we obtain 

\{u„,[d^,b{-e)]{b{0)a{ep){cm\li X {eH))))\ < C{Q){1 + S). 

To bound the last term on the right-hand side of (3.27), we use (2.6) to deduce 
that 

{u^,[d",b{-d)]{b{9)a{ep)div^f{eu))) 

<Co||u„|U2||[a",6(-0)](&(0)a(ep)divM/(eu))|U2 

<C{Q){\\b{-e)\\w,^^\\b{0)a{ep)div^{eu)\\H^ 

+ \\bi-9)\\H^^4b{9)a{ep)dw^{eu)U^) 
<C{Q){1+S). 

Hence, it holds that 

\{nc.,g2)\<C{Q){l+S). (3.28) 

Since 173 is similar to gi in structure, we easily get 

\{d.,93)\<C{Q){l+S). (3.29) 

Therefore, it follows from (3.26), (3.28)-(3.29), the positivity of &(-6'), and the 
definition of O, Oq, Q and S, that there exists a constant li > 0, such that for 
t e [0, T] and T = min{ri, 1}, 

sup II (p 

)('^)lll/= + i + '1 / ||V(Ua,^a)||_ff-s+i(T)dT 

Te[o,t] J 

<CiOo) + C{0{t))t + C{0{t)) f S{T)dT 

Jo 

<CiOo) + C{0{t))Vt 

< C{Oo)exp{VTCiOiT)}. 

Summing up the above estimates for all a with < |a| < s + 1 leads to the 
desired inequality (3.12). □ 

In a way similar to the proof of Lemma 3.3, we can show that 

Lemma 3.4. Let s > 4 and (p, u,H,6') be a solution to (1.12)^(1.15), (1.21) on 
[0,Ti]. Then there exists a constant I2 > 0, such that for any e G (0,1] and 
t € [0, T], T = min{ri, 1}, it holds that 

sup \\{e\,e^u,ei9-9)iT)\\H.^2+l2{ T ||V(e'u, e(?)||^.+.(T)dT}'^' 

<C(Oo)cxp{(VT)C(0(T))}, (3.30) 

where 0{T)) and Oq are defined by (2.8) and (2.9), respectively. 

Recalling Lemma 2.2 and the definitions of Q and S, one finds that 

||at(ep, eu,eH, 0)11^.-1 <C(Q), ||9*(ep, eu, eH, ^)|U. < C(Q)(1 + 5), (3.31) 
e\\dt{ep,eu,eli,9)\\H^ < C(S). (3.32) 
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Moreover, it follows easily from Lemmas 3.1-3.4 and the equation (1.15) that for 
some constant I3 > 0, one has 

sup \\edt0\\ls + h f \\Viiedt)efH.{T)dT<C{Oo)exp{VTC{0{T)}. (3.33) 

re[0,T] Jo 

3.3. H''~^ estimates on (divu,Vp). To establish the estimates for p and the 
acoustic part of u, we first control the term {edt){p, u). To this end, we start with 
a i^-estimatc for the linearized system. For a given state (poj Ug, Hq, ^o); consider 
the following linearized system of (1.12)-(1.15): 

dtp+ (uq • V)p + idiv(2u — Ra{epo)h{df))\/ 9) = ea(epo)[f' curlHo : curlH] 

+ ea(epo)«'(uo) : Vu + fia{epo)b{eo)Vpo ■ Ve + A, (3.34) 

b{-9o)[dtU + (uo • V)u] + -^1 = a(epo)[(curlHo) x H + div^(u)] + h, (3.35) 

e 

dtU - curl (uo X H) - PAH = /a, divH = 0, (3.36) 

dtO + (uo • V)6' + divu = e^a(epo)[i^curlHo : curlH + *(uo) : Vu] 

+ Ka(epo)div(6(0o)V(?) + /4, (3.37) 

where we have added the source terms fi (1 < i < 4) on the right-hands sides of 
(3.34)-(3.37) for latter use, and used the following notations: 

a{epo) := e~'P° , b{0o) -.^ e'" . 

The system (3.34)~(3.37) is supplemented with initial data 

(p,u,H,0)|t=o = (Pin(a;),Ui„(x),Hin(a;),0i„(a;)), x (3.38) 

Lemma 3.5. Let (p, u, H, 0) be a solution to the Cauchy problem (3.34)"(3.38) on 
[0,r]. Then there exist a constant I4 > and an increasing positive function C(-), 
such that 

sup ||(p,u,H)(T)||i3+/4 r||V(u,H)||i.(r)dr 

Te[o,T] Jo 

<e^^(^")|l(^,,u,H)(0)||i.+C7(i?o)e^^(««' sup ||V%)||i. 



Te[o,T] 



+ C(i?o) / ||V(6U,6H)||2,dr + C(i?o) / \\Ve\\l,{T)d7 
Jo Jo 

+ C{Ro) r {WhWh + \\f2\\h + \\f3\\h + llVMIi.} (r)dr, (3.39) 



whe 



i?o= sup {||9A(r)||L=c,,||(po,uo,Ho,0o)(r)|Ui.oo}. (3.40) 

-re[0,T] 



Proof. Set 

(p, u, H, 9) = {p, 2u - fia{epo)b{eo)Ve, H, 6). (3.41) 
Then p and H satisfy 

dtp + (uo • V)p + idivu = ea(epo)[i^curlHo : curlH] + ^a{epo)^{uo) : Vu 

6 ^ 
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+ |a(epo)*(uo) : V(Sa(epo)&(eo)V^) + Sa(epo)K^o)Vpo • ^0 + h (3.42) 

and 

dtii - curl (uo X H) - f^AH = /g, divH = 0, (3.43) 

respectively. One can derive the equation for u by applying the operator V to 
(3.37) to obtain 

dtV9 + (uo • V)V6' + ivdivu + ivdiv(Ka(epo)&(6'o)V6i) 

= V |e^a(epo)[;^curlHo : curlH]| + {fP' a{epo)'^ {vlo) : Vu} 

+ {e2a(epo)*(uo) : V(Sa(epo)6(0o)V^~)} 

+ V {Ka(epo)div(6(6lo)V^)} + [V, Uq] • V0 + V/4. (3.44) 
If we multiply (3.44) with iKa(epo)i we get 
\h{-eo) [dt{-na{epQ)h{eo)V~e) + (uo • V)[Sa(epo)6(eo)V0]} 
^ ^b{-e^)dt{a{epo)h{e^)}V~e + |6(-0o) {uo • V[a{epa)b{e^)]Ve] 
— — Ka(epo)Vdivu — — Ka(epo)Vdiv(Ka(epo)^(^o)V6') 

+ ^Ka(epo)V |e^a(epo)['^curlHo : curlH]| + i«;a(epo)V {e^a(epo)*(uo) : Vu} 
+ isa(epo)V |e2a(epo)*(uo) : V(Ka(epo)&(6'o)V^)} 
+ isa(epo)V {Ka(epo)div(6(6lo)V6')} 

+ ^Sa(epo)[V,Uo] • iKa(epo)V/4. (3.45) 

Subtracting (3.45) from (3.35) yields 
i5(-0o)[5tii + uo-Vii] + ^ 

= - '^b{-eo)dt{a{epa)h{eo)]Ve - ^^-^o) {uq • V[a(epo)6(^o)] V^} 
+ -Ka(epo)Vdivu + -Ka(epo)Vdiv(Ka(epo)&(^^o)V6') 

— -Ka(epo)V |e^a(epo)[f'CurlHo : curlH]| 

- hia{epo)\/ {e^a{epa)-^{no) : Vu} 

- iKa(epo)V {e2a(epo)*(uo) : V(Ka(epo)6(^o)V^)} 

— -Ka(epo)V |Ka(epo)div(fe(6'o)V6')| + a(epo)[f'CurlHo : curlH] 
+ ia(epo)[^'(uo) : Vu] + ia(epo)[*(uo) : V(Sa(epo)&(^o)V^)] 
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- iKa(epo)[V,Uo] • + a(epo)[(curlHo) x H] 

+ ia(epo)div*(Ka(epo)&(6'o)V^) + ia(epo)div*(u) - ^fia{epo)V fi + /a 

14 1 1 

■=Y.^^ + ;Ta(ePo)div*(u) - -fia{epo)V U + /a- (3.46) 



Multiplying (3.42) by p, (3.43) by H, and (3.46) by u, respectively, integrating the 
results over R"^, and summing them together, we deduce that 

d (1 

1 ,„ . , „ , . 1 



:{liP^P) + \{b{-Oo)n,h) + i(H, H)} + VH||i. 



((uo • V)p,p) + -{dtbi-eo)u,u) - -(5(-0o)(uo • V)u,u) 
{ea{epo)[D cmlHo : curlH],p) + |(a(epo)*(uo) : Vu,p) 
|(a2(epo)6(eo)*(uo) : V(V^),j5) + {Ka{epo)b{9a)V Po ■ VO,p) 



14 ^ 

^ fi) + -(a(epo)div*(u), u) 



2 

1=1 

- i (Sa(epo)V/4,u) + (/2,u) + (/3,H) + (A,p), (3.47) 

where the singular terms have been canceled out. 

Now, the terms on the right-hand side of (3.47) can be estimated as follows. 
First, it follows from the regularity of (po, uq, Hq, 6*0), a partial integration and 
Cauchy-Schwarz's inequality that 

^ \{dtb{~9o)u,u)\ < l\\dtb{-9o)\\L^\\u\\h < C{Ro)\\u\\h, (3.48) 



\{{uo-V)p,p)\^l 



(divuo)|p| dx 



<CiRo)\\p\\i., (3.49) 



i|(fe(-^^o)(uo•V)u,u)| <C(i?o)j|u||i2, (3.50) 
\{ea{epo)[9cmino : curlH],p)| < C{Ro){\\eVil\\l, + \\p\\l,), (3.51) 
||(a(epo)*(uo) : V{i,p)| < CiRo){\\eWu\\l. + \\p\\h), (3.52) 



--\{a^{epo)b{0o)^iu„) : V(V0),p)| < CiRo)\\p\\h 



+Gi(epo,eo) E ll<9"M)lii^' (3-53) 

|a|=2 

|(Sa(epo)fo(^?o)Vpo • V9,p)\ < C(i?o)(|| V^li. + ||p||i.), (3-54) 

where Gi(-, •) is a smooth function. Similarly, one can bound the terms involving 
hi as follows. 



J2m,u)\ <|||VH||i, + f llVulli, + f lldivulli, 

2—1 

+ C(i?o)|lii|li,+C(i?o)||V(eii,0~)||i^+G2(epo, 00)11 A^li^, (3.55) 
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where G'2(-, •) is a smooth function. 

For the dissipative term ^(a(epo)div^(u), u), wc can employ argmncnts similar 
to those used in the estimate of the slow motion in (3.22)-(3.24) to obtain that 

-l(a(epo)divvI/(u),u)>^(||Vu|li. + ||divu||i.)-C^(i?o)|jfl|li.. (3.56) 

Finally, it follows from (3.47)-(3.56), and Cauchy-Schwarz's and Gronwall's in- 
equalities that (3.39) holds. □ 

In the next lemma we utilize Lemma 3.5 to estimate ((e9t)p, (e9t)u, (e9i)H). 

Lemma 3.6. Let s > 4 and {p, u, H, 9) be the solution to the Cauchy problem 
(1.12)-(1.15), (1.21) on [0,r]. Set 

(p^,u^,H^,0^) := dP{{edt)p, {edt)u, {edt)H, {edt)0), 

where 1 < |/3| < s — 1. Then there exist a constant > and an increasing positive 
function C{-), such that 

sup ||(p^,u^,H^)(T)||i. +^5 / ||V(u^,H/?)|li2(r)dr 

re[0,T] Jo 

<C(Oo)exp{(%/T)C(0(T))}. (3.57) 
Proof. An application of the operator d^{edt) to the system (1.12)-(1.15) leads to 

dtPi3 + (u • V)p^ + -div(2u^ — Ra{ep)b{6)V 6 p) = ea(ep)[t?curlH : curlH^] 

+ ea(ep)^'(u) : Vu/j + Ra{ep)b{9)Vp ■ VOp + .gi, (3.58) 

b{~6)[dtnp + (u- V)u^] + ^ = a(ep)[(curlH) x + div*(ua)] + 52, (3.59) 

5(H/3 - curl (u X H;3) - i^AH^ = 53, divH/? = 0, (3.60) 
dtOp + (u • \/)e0 + divu^ = e^a(ep)[z?curlH : curlH^] 

+ e'^a{ep)'^{u) : Vug + Ra{ep)div(b(e)\/ef)) + 54, (3.61) 

where 

51 ~~[dP{edt),u]-\/p+^[d^{edt),{Ra{ep)b{0))]Ae 

+ i[9''(e5t), V(Sa(ep)6(6l))] - S/O + eD [d^ (edt) , a{ep)cuT\U] : curlH 
e 

+ e[df^{edt),a{ep)^{u)] : Vu + [d^{edt), Ra{ep)b{e)V p] ■ VO, (3.62) 
~g, :^ - [d^iedt),bi-e)]dtu ~ [df'iedt),bi-6)u] ■ Vu, 

+ [d^{edt),a{ep)cm\U] x H + [5'^(e9t), a(ep)]div*(u), (3.63) 

53 ■.=d^{edt){cm\{u x H)) - curl(u x H^), (3.64) 

54 ■.^-[d^^{edt),u]-We + e^D[d^^{edt),a{ep)cuT\lI] :curlH 

+ e^[d^iedt),a{ep)^{u)] : Vu 

+ Rd^^{edt){a{ep)div{b{9)\/9o,)) - Ka(e)div(&(6')V6'„). (3.65) 
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It follows from the linear estimate (3.39) that for some h > 0, 

sup ||(p0,U/3,H^)(r)!li.+/4 r!lV(u^,H^)|li.(r)dr 

rG[0,T] Jo 

<e^^(«)||(p^,u^,H^)(0)||i.+C(i?)e^^(«) sup ||V0^(r)||i. 

Te[o,T] 

+ TC{R) sup \\V{eu^,eHf;){T)\\l,+C{R) r \\Ve^\\j,,{T)dT 

re[0,T] Jo 

+ C{R) [ {||gi!li. + il,g2||i2 + i|,g3|li2 + i|V.g4||i2}(r)dr, (3.66) 



where R is defined as Rq in (3.40) with {po, Uq, Hg, ^q) replaced with {p, u, H, 

\l2, IMh, IMh, and |lVg4|li 



It remains to control the terms ||.gi||^2, 11.9211^2, 115311^2, and ||Vg4jl^2. The first 



term of gi can be bounded as follows. 

\\[d\edt),u]-Vp\\^, <eCo{MH^^4iedt)Vp\\H^^2 + \\iedt)u\\H^-4Vp\\Hs-^ 

<C{Q). (3.67) 
Similarly, the second term of gi admits the following bound; 

-\\[d^{edtl{Ra{ep)b{e))]A9\\ 
e 

< C4\\a{ep)b{e)\\H^-4dtAe\\Hs-2 + \\dt{a{ep)b{e))\\Hs-i\\Ae\\H^-.) 
<C{Q){l + S). (3.68) 

The other four terms in gi can be treated similarly and hence can be bounded by 
C(Q)(l + 5). 

For the first term of 52, one has by the equation (1.13) that 
[d^{edt),b{-e)]dtu =[d^{edt),b{-e)]{{u ■ V)u} 

+ -[d^iedt),b{-e)]{b-\-eo)Vp} 
e 

- [d^{edt),b{-e)]{b-\-9)aiep)[{cuT\lio) x H]} 

- [d^iedt),bi-e)]{b-\-9)a{ep)div^{u)}. (3.69) 

Note that the terms on the right-hand side of (3.69) have similar structure as that 
of gi. Thus, we see that 

\\[df^{edt),b{~e)]dtu\\L2 < C{Q){1+S). (3.70) 

Similarly, the other four terms of g2 can be bounded by C(Q)(1 + 5). 
Next, by the identity (2.4), one can rewrite §3 as 
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(eaO,divu]H- [a^(e5t),u] • VH + ^[9^(eat),Vu,]H. (3.71) 



Following a process similar to that in the estimate of gi , one gets 

\\h\\L-2<Cmi+S). (3.72) 

And analogously, 

llV.g4!lL2 <C(Q)(l+5). (3.73) 
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We proceed to control the other terms on the right-hand side of (3.66). It follows 
from (3.33) that 

C(i?)e^'=^(^) sup \\S/9p{T)\\l, <C{0)exp{VTC{0{T)} 
Te[o,T] 

and 

\\AOp\\l^_iT)dT < [ \\ie^t)9\\l.^^iT)dr < CiOo)exp{VfCiOiT))}. 



Thanks to (3.32), one has 

TC{R) sup ||V(eu0,eH^)(T)||i. <TC(0(r)) sup \\{edt){eu,eH){r)\\%. 
Te[o,T] Te[o,T] 

< TC{0{T)). (3.74) 

Then, the desired inequality (3.57) follows from the above estimates and the in- 
equality (3.66). □ 

Now we are in a position to estimate the Sobolev norm of (divu, Vp) based on 
Lemma 3.6. 

Lemma 3.7. Let s > 4 and (p, u, H, 6) he the solution to the Cauchy problem 
(1.12) -(1.15), (1.21) on [0,r]. Then there exist a constant /g and an increasing 
function C{-), such that 

sup {||p(t)||h., + ||divu(T)|lH.-i} + /6 / {llVpIl^. + |lVdivu|l^.-i}(T)dT 
Te[o,T] Jo 

< C(Oo)exp{(VT + e)C(0(r))}. (3.75) 
Proof. Rewrite the equations (1.12) and (1.13) as 

divu = - ]^{edt)p - |(u • V)p + idiv(Sa(ep)6(0)V0) + ^a(epo)|curlH|2 

+ — a(ep)«'(u) : Vu + —a[ep)b{e)Vp ■ V6', (3.76) 

Vp = - b{-e){edt)u - eb{-e){u ■ V)u 

+ ea(ep)[(curlH) x H] + ea(ep)div*(u). (3.77) 

Then, 

||divu|U.-i <Co\\{e^t)p\\H'-^^Cae\\xl\\Hs-^\\Vp\\Hs-. 

+ Co||div(Ka(ep)&(0)V^?)||ffs-i +Co||a(epo)|lL-!lecurlHjl2,_i 
+ Co||a(ep)||Loc|j*(eu) : (eVu)||^=-i 

+ C4a{ep)b{0)\\LA{^Vp)\\Hs-A\V9\\Hs^.. (3.78) 
It follows from Lemmas 3.2-3.4 and 3.6, and the inequalities (3.31)-(3.33) that 
\\{edt)p\\H^-i <C(Oo)exp{(VT-t-e)C(0(T))}, 
e||u||^.-i||Vp|U.-i <eC(0), 
\\dw{Ra{ep)b{e)Ve)\\H.-. <Co|| A(?||^.-i + Co|| V0||h=-i 

<C7(Oo)exp{(x/T + e)C(0(T))}, 
||ecurlH|l^.-i <C(Oo) exp{(VT + e)C(0(T))}, 
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||*(eu) : (eVu)|U..-i <C(Oo) exp{(VT + e)C(0(T))}, 

||(eVp)||H-^-i||V0||H^.-i <C{Oo)exp{{VT + e)C{0{T))}. 

These bounds together with (3.78) imply that 

sup \\divu\\{T)\\Hs-i<C{Oo)exp{iVT + e)C{0{T))}. (3.79) 
Te[o,T] 

Similar arguments applying to the equation (3.77) for Vp yield 

sup \\p\\H~^+k[ \\\/ p\\ j,s {T)dT <C{Oo)cxp{{VT + e)CiOiT))} (3.80) 
Te[o,T] Jo 

for some positive constant /g > 0. 

To obtain the desired inequality (3.75), we shall establish the following estimate 

2 



||Vdivu||^.-i(r)dr < C(C'o)cxp{(VT + e)C(C'(T))}. (3.81) 

To this end, for any multi-index a satisfying 1 < |q;| < s, one can apply to 
(3.75) and then take the inner product with 3"divu to obtain 

r |ja"divu||i. (r)dr = - i r {d"{edt)p,d''dwu){T)dT 

+ / (S,a"divu)(T)dT, 

Jo 



(3.82) 



where 



S = - |(u • V)p + ^div(Sa(ep)6(0)V0) + ^a(epo)|curlH|2 

2 - 

+ y a(ep)^'(u) : Vu + ^a{ep)b{e)Vp ■ V9. 

It thus follows from (3.57) and similar arguments to those for (3.80), that for all 
1 < |a| < s, 

[ \\d^Eir)\\l,dT < C{Oo) exp{{VT + e)C{0{T))}, 
Jo 

TkS, a" divu) I (r)dr 
Jo 



whence, 



<C(Oo)cxp{(VT + e)C(0(T))}<^ / ||9"divu||i. (r)dr 



T 

2 



1/2 



< C(Oo)exp{(VT + e)C(0(T))} + i^ |1 9"divu|l (T)dT. 

For the first term on the right-hand side of (3.82), one gets by integration by 
parts that 

cT IT 



1 I (a"(e(9f)p,<9"divu)(T)dr = - e9"div(eu)) 

2 Jo 2 

(9"Vp,9"(eat)u)(T)dT. 
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By virtue of the estimate (3.12) on {eq, eu, — 9) and (3.80), we find tliat 



i(a>, ea"div(eu))'^ 



(9"Vp,a"(eat)u)(T)dr 



< sup {|lp|lHHleu||ff=+i} 

re[O.T] 



<C7(Oo)exp{(VT + e)C(0(T))}, 



/ WVp\\i.{T)dT 



1/2 



^ W{edtMh{T)dT' ' 



<C(Oo) cxp{(VT + e)C(0(T))}. 

Tfiese bounds, together with (3.82), yield the desired estimate (3.81). This com- 
pletes the proof. □ 

3.4. H^^^ estimate on curlu. The another key point to obtain a uniform bound 
for u is tlic following estimate on curl u. 

Lemma 3.8. Let s > 5/2 be an integer and (p, u,H,0) be the solution to the 
Cauchy problem (1.12)-(1.15), (1.21) on [0,T]. Then there exist a constant Iq > 
and an increasing function C(-), such that 

sup {||curl(6(-0)u)(r)||^._, +||curlH(r)||^._,} 
Te[o,T] 

+ k r {||Vcurl(&(-0)u)ll^._i + !lVcurlH(T)||^._i }(T)dT 
Jo 

<C{Oo)cxp{{VT + e)C{0{T))}. (3.83) 

Proof. Applying curl to the equations (1.13) and (1-14), using the identities (2.1) 
and (2.2), and the fact that curlV — 0, one infers 

at(curl {b{-d)u)) + (u • V)(curl {b{-9)u)) 

= curl {a(ep) (curl H) x H} + /2div{a(ep)&(6l)V(curl (6(-6')u))} + Ti, (3.84) 

5t(curlH) - curl [curl(u x H)] - i>A(curlH) = 0, (3.85) 

where Ti is defined by 

Ti --fldiY {a{ep){Vb{9)) curl {b{-0)u)) - /2Va(ep) • V(6(6i)curl (6(-e')u)) 

- fia{ep)A{(yb{9)) x (5(-6')u)) - Va(ep) x (^Au + [fl + A)Vdivu) 

+ curl {b{-9)udt9) + [curl , u] • V(6(-6i)u) + curl (b(-9)u{u ■ WO)). 

For any multi-index a satisfying < |q;| < s ~ 1, we apply 9" to (3.84) and 
(3.85) to obtain 

9(9" (curl (6(-6l)u)) + (u • V)[a"(curl (6(-6l)u))] 

= 9"curl{a(ep)(curlH) x H} 

+ fLdiv{a{ep)b{9)W[d"'{cm\{b{-0)u))]} + d^'Ti + T2, (3.86) 

9(9" (curl H) - 5" curl [curl (u x H)] - i?A(curlH) = 0, (3.87) 

where 

T2 := - [9", u] • V[5"(curl {b{-9)u))] 
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- [9",div(a(ep)fe(6'))]V[a"(curl(fe(-6')u))]. 

Multiplying (3.86) by d°'{cm\ {b{—9)u)) and (3.87) by c)"(curlH) respectively, sum- 
ming up, and integrating over R'^, we deduce that 

i|{||5"(curl(fo(-^^)u))||i. + i|9"(curlH)||i.} + ^^||a"(curlH)||i. 

+ (a(ep)6(6l)V [(9" (curl (fe(-6l)u))], V[5"(curl (6(-6')u))]) 
= - ((u • V)[a"(curl (&(-6l)u))], a"(curl (6(-6')u))) 

- (9"curl{a(ep)(curlH) x H}, 0"(curl (fo(-6')u))) 

+ (9"curl [curl(u x H)], ^"(curlH)) + (9"Ti + T2, ^"(curl (&(-6')u))) 

:= Ji +J2 + J3 + J4, (3.88) 

where J7i (i = 1, • • • ,4) will be bounded as follows. 
An integration by parts leads to 

\Ji\ < ||divu|Uoo||V[a"(curl(&(-0)u))]||i.. 

By virtue of (2.3), the Cauchy-Schwarz's inequality and Moser-type inequality (see 
[33]), the term J2 can be bounded by 

\J2\ <|(a"{a(ep)(curlH) x H}, 9"curl (curl (6(-6')u))) | 

<|a"{a(ep)(curlH) x H}||i2||a"V(curl (fe(-6')u))||i2 

<m\\d''V{cm\{b{-e)u))\\l. 

+ C„{||curlH||i^||a(ep)H||^._, + ||a(ep)H||i^||curlH||^._,}, 

where 772 > is a sufficiently small constant independent of e. 

If we integrate by parts, make use of (2.3), curl curl a = Vdiva — Aa and the 
fact that divH = 0, we see that the term J/3 can be rewritten as 

J3 = (a"curl (u X H), 9" AH) , 

which, together with the Moser-type inequality, implies that 

m<CiS)+V3\\il^{r)\\l^„ 

where 7^3 > is a sufficiently small constant independent of e. 

To handle J74, we note that the leading order terms in Ti arc of third-order in 9 
and of second-order in u, and the leading order terms in T2 arc of order s + 1 in u 
and of order s -I- 1 in (ep, 9). Then it follows that 

IJ4I < Co(||9"TilU. + l|T2|U.)l|a"(curl(6(~0)u))|U. 
< C(5)|ia"(curl(6(-0)u))|li.. 

Putting the above estimates into the (3.88), choosing 772 and 773 sufficient small, 
summing over a for < |a| < s — 1, and then integrating the result on [0,T], we 
conclude 

sup {||curl(6(-^)u)(r)||^.„, -M|curlH(r)||^._,} 

-rG[0,T] 

+ h I {llVcurl (6(-0)u)|l^._i + |lVcurlH(T)|l^._i }(T)dr 
Jo 

< Co{\\cm\{h{-9)n){Q)\\H^-^ + i|curlH(0)||l,.„i} 
+ C{Oo) exp{( + e)C{0{T))} 
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< C{Oo) exp{(VT + e)C{0(T))}. 

□ 

Proof of Proposition 2.4- By the definition of the norm j| • \\s.r, and the fact 

||v||^™+i < A'(||divv||H„. + j|curlv||ff™ + ||v||ff..), V V e H'^+^R^), 
Proposition 2.4 foUows directly from Lemmas 3.1, 3.3, 3.7 and 3.8. □ 

Once Proposition 2.4 is established, the existence part of Theorem 1.1 can be 
proved by applying the arguments in [2,39] directly, and hence we omit the details 
here. 

4. Decay of the local energy and zero Mach number limit 

In this section, we shall prove the convergence part of Theorem 1.1 by modifying 
the arguments developed by Mctivier and Schochet [39], see also some extensions 
in [1,2,37]. 

Proof of the convergence part of Theorem 1.1. The uniform estimate (1.24) implies 
that 

sup |l(p^u^H^)(r)||ff. + sup ||^?^ - ^||^e+i< +00. 

Te[0,To] re[0,To] 

Thus, after extracting a subsequence, one has 

(p%u^)-(p,w) weakly-* in L°°(0, Tq; i/^CM^)), (4.1) 

H'^B weakly-* in L°°{Q,To]W{R^)), (4.2) 

e^^e^d-e weakly-* in L°°{Q,Tq]H''+'^{R?)). (4.3) 

It follows from the equations for and 6'^ that 

9tH^c)tr eC([0,To],i?-'-2(R3))^ 

Hence, after further extracting a subsequence, we obtain that for all s' < s, 

W^B strongly in C{[0,To], hC{^^)), (4.4) 

B'-B^d-e strongly in C{\0,TqI H^^+\R?)), (4.5) 

where the limit B e C{[0,To\, H^^^{M?)) n L°°(0, Tq; i7f„^(R3)) and {'d - B) e 
C([0,To],if,;'+HM3)) n (0, To; i7f 
Similarly, from (3.83) we get 

curl (g-'^'u') ^ curl (e-'^w) strongly in C([0, Tq], H^~^ [R^)) (4.6) 

for all s' < s. 

In order to obtain the limit system, one needs to show that the limits in (4.1) 
hold in the strong topology of ^^(O, Tq; i/fJj.(R^)) for aU s' < s. To this end, we 
first show that p = and div(2w — Re^V'd) = 0. In fact, the equations (1.12) and 
(1.13) can be rewritten as 

e dtp' + div(2u'^ - Ke-'P'+'^'\/B') = e/% (4.7) 
ee-^'^tu' + Vp'^eg". (4.8) 
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By virtue of (1.24), f and are uniformly bounded in C([0, To], H^-^M.^)). Pass- 
ing to the weak limit in (4.7) and (4.8), respectively, we see that Vp = and 
div(2vir - Se^'Vi?) = 0. Since L°°(0, Tq; iJ^(M3)), we infer that p = 0. 

Notice that by virtue of (4.6), the strong compactness for the incompressible 
component of u*^ holds. So, it is sufficient to prove the following proposition 
on the acoustic components in order to get the strong convergence of u*^. 

Proposition 4.1. Suppose that the assumptions in Theorem 1.1 hold. Then, p"^ 
converges toQ strongly in L^{0,To; Hf^^{M.^)) and div(2u'^— Ke~'P VS'^) converges 
to strongly in L'^{0,Tq; H^'^-'^{W)) for all s' < s. 

The proof of Proposition 4.1 is based on the following dispersive estimates on 
the wave equation obtained by Metivier and Schochct [39] and reformulated in [2]. 

Lemma 4.2. ( [2,39]) LetT >Q andv' be a bounded sequence inC{[Q,T],H'^{M?)), 
such that 

e^dt{a'dtv')^V ■ ib'Vv') = c', 

where c'^ converges to strongly in L^(0, T; L^(R'^)). Assume further that for some 
s > 3/2 + 1, the coefficients (a'^,6^) are uniformly bounded in C([0, T], iJ*(R'^)) and 
converge in C([0, T], i7['^^(R'^)) to a limit (a, 6) satisfying the decay estimate 

\a{x,t)-a\<Co\x\-'-'^, |V,a(x, t)| < Co|x|-2-C, 

\b{x,t) - 6| < Co|x|-i-«, |V,6(x,t)| < Co|x|-2-C, 

for some positive constants a, b, Cq and ^. Then the sequence v'^ converges to 
strongly m L^{0,T; Lf^^{R^)). 

Proof of Proposition 4-1- We fist show that converges to strongly in L^{Q, Tq; 
Hiici^^)) for aU s' < s. Applying e'^dt to (1.12), we find that 

e^dtidtp" + (u^ • V)p"} + e9t{div(2u^ - Re-^P'+^'VO^)} 

= t'dtie^'P^ [z>|curlH'|2 + *(u^) : Vu']} + e^dtihe^'^^+^'Vp' ■ VO'}. (4.9) 

Dividing (1.13) by and then applying the operator div to the resulting equa- 
tion, one gets 

e^tdivu" + div(e^'' Vp<^) = - ediv{(u'= • V)u'} 

+ ediv{e-"P'+^'[(curlH) x H + div5''=(u'^)]}, (4.10) 

Subtracting (4.10) from (4.9), we have 

-div(e^Vp^) =eF^(p^u^H^0^), (4.11) 

where i^'^(p', u*^, H*^, 6"^) is a smooth function in its variables with F(0) = 0. By the 
uniform boundcdness of (p*^, u*^, H*^, O"^) and the strong convergence of 0*^, one infers 
that 

F'{p',u\-H.\e')^Q strongly in L^{Q,To]L^{M.^)), 

and the coefficients in (4.11) satisfy the conditions in Lemma 4.2. Therefore, we 
can apply Lemma 4.2 to obtain 

/^O strongly in L^^q, Tq; iL(R'))- 



LOW MACH NUMBER LIMIT FOR THE FULL COMPRESSIBLE MHD EQUATIONS 25 



Since p'^ is bounded uniformly in C {[0,To], H'^ (R^)) , an interpolation argument 
gives 

strongly in 2.^(0, To; iJfo'^(M^)) for all s' < s. 

Similarly, we can obtain the strong convergence of div(2u'^ — k'^c"'^ VO'^). This 
completes the proof. □ 

We continue our proof of Theorem 1.1. It follows from Proposition 4.1 and (4.5) 
that 

divu"->divw strongly in L^{0,Tq;H^^^'^{W^)). 

Thus, using (4.6), one obtains 

^ w strongly in L^iQ, Tq; i7i;'^(M^)) for aU s' < s. 

By (4.4), (4.5), and Proposition 4.1, it holds that 

Vu^-^Vw strongly in L^{0,To; H^^;:^{R^)); 

VB strongly in L^{0, Tq; H^'^-\R^)); 
V9' strongly in L^{0, Tq; H^^^r^R^)). 

Passing to the limits in the equations for p^, H'^, and 6*^, we see that the limit 



(0,w, B,i?) satisfies, in the sense of distributions, that 

div(2w - K = 0, (4.12) 

dtB - cmi (w X B) - PAB = 0, divB = 0, (4.13) 

dt^ + (w • V)d + divw = fidiv{e^Vi9). (4.14) 



On the other hand, applying curl to the momentum equations (1.13), using 
the equations (1.12) and (1.15) on p'^ and 9'^, and then taking to the limit on the 
resulting equations, one sees that 

curl {dtie^'^w) + div(we"'' «) w) - (curlB) x B - div$(w)} = 

holds in the sense of distributions. Therefore it follows from (4.12)-(4.14) that 

e'^[dtw + (w • V)w] + Vtt = (curlB) x B + div$(w), (4.15) 

for some function n. 

Following the same arguments as those in the proof of Theorem 1.5 in [39], we 
conclude that (w,B,i?) satisfies the initial condition 

(w,B,^)|,=o = (wo,Bo,^9o). (4.16) 

Moreover, the standard iterative method shows that the system (4.12)~(4.15) with 
initial data (4.16) has a unique solution (w*,B*,i9*) e C([0, Tq], i?"(M''^)). Thus, 
the uniqueness of solutions to the limit system (4.12)-(4.15) implies that the above 
convergence holds for the full sequence of (p*^, u*^, H*^, 6''^). Therefore the proof is 
completed. □ 
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